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ABSTRACT 


In recent times, the most widely used trajectory modelling for artillery shells has 
been via mathematical models such as point mass model, modified point mass model and 
six-degrees-of-freedom model. Applications of these models require an a priori 
postulation of equations of motion governing the shell trajectory and to solve these 
equations, one needs reliable estimates of aerodynamic coefficients. Due to various 
assumptions used in arriving at the mathematical models and also due to non availability 
of rehable estimates of aerodynamic coefficients required to solve these equations, an 
alternative approach of using general function approximation capability of the feed 
forward neural networks (FFNNs) for estimating shell performance is explored. In the 
literature, no work has been reported wherein FFNNs have been used to develop neural 
models for shell performance. The present work addresses this aspect by way of 
proposing three distinct neural models for predicting trajectory variables. The models are 
validated for 155 mm Bofors shell HE77B data supplied in the form of range tables by 
AEDE, Pune. The estimated trajectory parameters like the range, the firing angle, the 
time of flight, etc., via the proposed neural models compare well with those listed in the 
supplied range tables. The neural models developed take into account variations in shell 
mass and its muzzle velocity, and variable atmospheric conditions like head/tail wind, 
crosswind, temperature and density that might prevail at the time of firing. It is shown 
that the proposed models can accurately predict i) the range obtainable for varying firing 
angles under prevailing atmospheric conditions, ii) the firing angle required to achieve 
desired range for known atmospheric conditions, and iii) the standard range that the shell 
would have achieved under standard atmospheric conditions. 
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CHAPTER 1 
INTRODUCTION 


Artillery forms an important wing of the army to provide firepower during war as 
well as during cross-boarder skirmishes with the enemy. Artillery generally falls into 
three basic categories: guns, howitzers and mortars. The principle difference between 
them being the trajectory of the round fired. A gun has a high muzzle velocity and a very 
flat trajectory. Normally a gun is used in a direct fire mode where the target can be seen 
and penetration is desirable. Howitzers have a somewhat lower muzzle velocity and arc 
their shells onto a target. They are used in both a direct fire and indirect fire mode. This is 
especially useful when an enemy is concealed behind a prepared position or the artillery 
men desire to have a shell explode over an enemy’s head. The air-burst does less damage 
to hardened targets, but causes many more human casualties due to shell fragmentation 
covering a large area. Mortars have a very pronounced arc of flight. They have a 
relatively low muzzle velocity and are unsuitable for direct fire. Their principle value 
comes from being able to lob shells behind an obstacle, such as a fortification or a hill. 
They are not very accurate and dependent upon the amount of propelliag powder to 
determine the point of impact. 

The effectiveness of artillery is largely judged by the accuracy m hitting the 
targets. The accuracy and reliability is influenced by the design criteria used in designing 
the shell. Artillery shells are a class of projectiles around which much of the aerobaUistic 
theory was origmally developed, and they continue to form a significant part of the 
aeroballistician's interest. 
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Study of the motion of projectiles through an external medium is known by the 
name of external ballistics. If the external medium is the earth's atmosphere, then external 
ballistics become synonymous with aeroballistics. 

Each specific projectile will have its own individual design criteria with regard to 
payload, range, accuracy, etc. The practical design criteria applying to a conventional 
artillery shell which most affect the aeroballistic aspects are: 

1) High payload carrying potential 

2) Range at high charges (muzzle velocity around Mach 2-3) must be maximized - 
hence low drag 

3) Accuracy and repeatability 

4) Ease of manufecturing. 

The accuracy of artillery shell depends on the following major factors: 

1) Muzzle velocity irregularity due to variation in charge 

2) Jump and throw-off: This occurs due to recoiling effect of gun. Due to this action, 
there is a vertical component, known as jump, and there is a horizontal component of 
that force known as throw-off. 

3) Meteorological effects like ambient temperature, density, head/tail wind, crosswind, 
etc. 

4) Differences between shells due to shape, size, mass, etc. 

5) Difference in yawing behavior, at the time of exit from the barrel, due to muzzle blast 
, and wind consideration. 
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The first two fectors are outside the control of the aeroballistician. The rest can be 


minimized by ensuring that the shell is both adequately gyroscopically and dynamically 
stable under all launch conditions, and has low drag. 

The conventional approach hitherto for understanding the in-flight behaviour of 
projectiles was to develop mathematical models that could predict all elements of the 
trajectory from launch to target. To this purpose, it becomes essential that all forces, 
moments and other terms (e.g., coriolis force due to the rotation of the earth) affecting the 
flight of the projectile are accounted for in a well defined mathematical form\ Beginning 
with the most simple but relatively inaccurate mathematical model, the in-vacuo 
trajectory model, more and more sophisticated models of increasing accuracy such as the 
point mass model, the modified point mass model and the sbc-degrees-of-freedom model 
have been developed. A brief description of these models is given in chapter 3. Presently, 
we only wish to point out that even the best of these proposed models have their 
limitations due to their inability to model all the problem variables adequately. For 
example, the initial conditions at the time of shell leaving the barrel are not accounted for 
by any of the proposed mathematical models. Furthermore, the models require 
aerodynamic coefiicients as input (e.g., drag coefficient, damping in roll derivative Qp, 
lift curve slope Clo, etc.) and the estimates available for these coefficients are not so 
reliable. Finally, the varying atmospheric conditions over the height traversed by the 
projectile are accounted for only in an adhoc manner by using one single weighted mean 
value of the variables like temperature, density, head/tail wind and crosswind. 
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It is thus realized that even the best of the mathematical models available to date 


are not reliable for field application because accurate predictions are difficult for; 1) the 
range obtainable for various firing angles, 2) the firing angle required for specified range. 

The limitations of the mathematical models so far used for predicting the 
performance of artillery shell motivated us to look at an alternative approach to 
modelling. The feed forward neural network^’^’"* provides one such potential way of 
modelling. 

The neural networks have been successfully used in such diverse fields as signal 
processing, pattern recognition, system identification and control. In recent years, neural 
models of aircraft aerod}mamics have been successfully developed for many applications 
relevant to Aerospace Engineering^’^’^. For example, neural modelling has been used for 
estimating aircraft stability and control derivatives of stable*’®’^”, unstable^ \ aeroelastic 
aircraft^^. It was envisaged that a neural model can be developed to replace the use of 
hitherto used mathematical models for solving the shell trajectory related problems. 

For real life situations, it was realized that the problem of shell trajectory, broadly 
speaking, could be divided into three categories: Under prevailing ambient atmospheric 
conditions 1) to predict range for the chosen fifing angle, 0, 2) to predict required firing 
angle to achieve specified range, 3) to predict range that would have been achievable 
under standard atmospheric conditions, using the range data obtained under existing 
atmospheric conditions. The next step is to identify the relevant set of input-output 
variables for the network for each of these problems. Finally, a suitable architecture for 
the Neural Networks is to be searched to achieve acceptable functional mapping between 
the input-output variables for each of the problem. 
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In the present work, all the above problems are addressed and adequately solved. 
For demonstrating the prediction capability of neural models developed, the data used is 
for the 155mm Bofor shell, supplied by ARDE, Pune. Details of the data supplied and its 
use for testing various neural models is given in chapter 3. The details of modelling for 
the three categories of problems related to the shell trajectory, along with results obtained 
from each of three models are discussed in chapter 4. A brief overview of the neural 
network precedes these chapters 3 and 4, and is given in chapter 2. The dissertation ends 
with chapter 5 containing conclusions and a few suggestions for the future work. 
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CHAPTER 2 

ARTIFICIAL NEURAL NETWORK 


2.1 Introduction 

A neural network is a parallel distributed processor that has a natural 
propensity for storing experimental knowledge and making it available for use. It 
resembles the brain in two respects, 1) Knowledge is acquired by the network 
through learning process, 2) Inter neuron connection strengths known as synaptic 
weights are used to store the knowledge. Neural networks are also referred to in 
the literature^’^’^'^ as neuro computers, connectionist networks, parallel distributed 
processors, etc. 

Artificial neural networks consist of group of neurons arranged in a 
layered structure. Each neuron receives signal firom the neurons in the layer 
previous to itself and passes a signal on to the neuron in the following layer. The 
relationship between the summed inputs to a neuron and its output is governed by 
an ‘activation function’^. Some of the commonly employed activation functions 
are the step fimctions, the tangent-hyperbolic function, and the logistic 
(sigmoidal) function. 

Out of these, the most oflen-used activation function is the sigmoidal 
function defined as 


F(x) = 


1 


1 + e 


-x/X 


where X is the logistic gain. 
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The sigmoidal function is continuous, monotonically increasing and continuously 
differentiable, and it asymptotically approaches fixed finite values as the input 
approaches infinity (+ or -). 

Amongst the artificial neural networks, the feed forward neural networks 
(FFNNs) have foimd the favour with most researchers for applications in 
aerospace engineering problems^’^’^. The feed forward neural network consists of 
source nodes that constitute the input layer and one or more hidden layers and an 
output layer. Feed forward neural networks have neurons arranged in layers like 
directed graphs, implying a unidirectional flow of signals, and thus are static m 
nature. This kind of modeling develops mput-output relationship of a black-box 
kind. Each of the connection between neurons is assigned its individual weight 
and it is adjusted so as to yield the required output corresponding to the known set 
of inputs. Assignment of weights is done during the training sessions of the 
network. 

2.2 Ba'ck Propagation Algorithm 

One of the efficient methods for training the FFNN is the back 
propagation algorithm (BP A). Back propagation algorithm consists of a forward 
and backward pass through different layers of the network. During the forward 
pass, the input vector is applied to the input nodes of the network and its effect 
propagates through the network layer by layer. The connective weights are all 
kept fixed during the forward pass. On the other hand, in the backward pass, the 
weights are updated in accordance with error correction rule. Specifically, the 
actual response of the network is subtracted from the desired response to produce 
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an error signal. This error signal is then propagated backward against the direction 
of the connective weights. There are the four steps for back propagation 
algorithm. 

a) Initialization: Start with a reasonable network configuration and set all the 
synaptic weights randomly. 

b) Forward Computation: Let a training example represented by x(n) be applied to 
the input nodes. Then the output of the network is computed by proceeding 
through the network, layer by layer. Next, the error signal is computed using 
difference between the desired response vector and the output vector. 

c) Backward Computation: The back propagation algorithm is based on 
optimizing a suitably defined error function. At each point, the local output error 
cost function defined by the sum of the squared error is computed. The weights of 
the network are adjusted in such a way that the mean squared error (MSE) is 
minimized. 

The MSE is given by 

1 m n _ 

MSE=— ^ Z Z[Yi(j)-Xi(j)p 

mxni^lj^l 

where Y and X are the desired and predicted outputs, n is the number of data 
points and m is the number of output variables. 

d) Step (b) to (c) are repeated for each training pair in the training set until the 
error for the entire set is less than the prescribed value or the number of iterations 
exceed the prescribed limit. 
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The training algorithm is recursive in nature and it needs repetitive 
training sessions to achieve the required learning. There are many network 
influential (tuning) parameters like the learning rate, the momentum rate, the 
number of hidden layers, the number of iterations, the logistic gain, etc., that 
affect the accuracy of functional mapping between the input and the output 
variables. There are no set rules for fixing values of these influential parameters. 
In literature, a few guidelines are available to guide the choice of these 
parameters. However, the final choices of fine tuning these parameters is to be 
achieved by trial and error for the given problem, and it is a crucial step in finding 
a suitable neural model for the problem Some of the guidelines and thumb rules 
for traming of neural network are given in Appendix B. 

The neural network model is first trained on the known sets of input- 
output pairs of e}q)erimental or recorded data. The trained network is then capable 
of predicting the required output based on known (measured) input variables. This 
approach does not require a mathematical model or a transfer function relating the 
input-output data but only a sufficient set of input-output pairs of data. The choice 
of inputs that are likely to affect the output to be predicted by the neural model is 
of critical importance. The selection of inputs for the given problem is based on 
the understanding of the physics of the problem and engineering judgement. 

In the present work, a neural model is proposed for the trajectory 
modelling of an artillery shell fired from a gun. A non-linear relationship is 
mapped by the neural network between the input variables such as range, muzzle 
velocity, wind velocity, ballistic air temperature (BAT), ballistic air density 
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(BAD) with the output variables such as time of flight (T), bearing corrections for 
drift (4^) and angle of firing (0). A schematic of such a neural model is shown in 
fig.l. For atypical neural model used in the present work, Fig.l, the choice of the 
input and output variables for mapping projectile dynamics is dictated by the 
physical understanding of the phenomenon governing the shell dynamics. Once 
the input and output variables for the feed forward artificial network (FFNN) are 
selected, the FFNN model of shell dynamics is achieved without the need of a 
formal model structure formulation. 


Input Hidden Output 

Layer Layer Layer 


Range 
Muzzle Velocity 
H/TWind 


Crosswind 

BAT 

BAD 

PW 



Fig. 1 Schematic of feed forward neural network. 
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CHAPTERS 


TRAJECTORY AND RANGE MODELLING OF ARTILLERY 

SHELL 

3.1 General 

Artillery shells are a class of projectiles around which much of the aeroballistic 
theory was originally developed. Even today, it still continues to be of interest and further 
investigations for many aeroballisticians and user agencies. In the present work, we are 
specifically to focus on flight of a spin stabilized, dynamically stable, conventional 
artillery projectile, possessing either rotational or some form of mirror symmetry. For 
such projectiles, there are two major factors that differentiate their aerodynamics from the 
classical aerodynamics of, say, aircraft aerodynamics: 

1) The geometric symmetry of the projectile implies that many aerodynamic terms are 
themselves symmetric. For example, for a shell with cruciform tails will have 
aerodynamic derivatives Cmia = -Cnp and Cniq = Cut. 

2) The second factor to be considered is due to spin rate imparted to many shells, giving 

rise to aerodynamic effects that are unique to aeroballistics. For example, the magnus 
forces and moments are hardly of concern for aircraft dynamics, but may play 
significant role for shell dynamics. 

Presently, the conventional approach to modelling of trajectory and range of 
artillery shell has been by postulating a suitable mathematical model consisting of 
equations of motion. There are numerous such forms of trajectory model involving 
different basic assumptions and having different conq>lexities of solution. We shall 
briefly outline, in an increasing order of complexities and accuracy, the more commonly 
utilized analytical trajectory models used at present time. 
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3.2 The in- vacuo trajectory model 

In the in-vacuo trajectory model all aerodynamic forces and moments acting on 
the projectile are neglected. It is thus highly inaccurate for all but the shortest-range 
trajectories and the lowest drag projectiles. However, it does have its uses, in defining 
concepts such as trajectory rigidity, and also for setting limits on parameters such as 
range and vertex height. The simplicity of the model renders analytical solution easily. 

3.3 The point mass model 

In this model, it is assumed that the only aerodynamic force acting on the 
projectile is drag. It provides fairly accurate estimates of range for adequately stable 
projectiles and can also be used to estimate the first order effects of wind. The point mass 
model is given as below. 


d^x 

dt2 


Ttpd^Cp 

8m 


dx 

dt 


■W. 


d^y 

-^ = -8- 
dt^ 


Tcpd^Cp J dy 


8m 


dt ^ 


d^z 


Ttpd^Cp 

8m 


fdz ) 

Vdt J 


(la) 


(lb) 


(Ic) 


where x denotes range, y denotes height, z denotes drift and Wx, Wy, Wzare, respectively, 
the X, y and z components of the wind velocity W. This model, however, does not 
account for the effect of spin of the shell and thus feils to predict the drift due to magnus 
moment. Also, it neglects the lift forces acting on the shell. 
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3.4 The Modified point mass model 

The modified point mass model is also known as 4 degree-of-freedom model (3 
spatial degrees-of-freedom plus axial spin). Its basis is a conventional point mass, in 
addition, the instantaneous equilibrium yaw is calculated at each time step along the 
trajectory so as to provide estimates of yaw, drag, drift and magnus force effects resulting 
firom the yaw of repose. The modified point mass equations of motion are given as 


dt 



ttr- 


Tipd^ 

16m 


%r,ry X v) - g^ — f + 2(o X u) 

PW' 


(2a) 


dt 16Ix P 


a r = — 


8p^x [vx (du/dt)] 
Ttpd^CMa 


(2b) 

(2c) 


The quantities in the above equation are defined as follows: 


dx 


u = — , v = u-W, r = x-R, R = (0,-R,0) [W denotes the wind vector] 

dt 


<a = (-Qcos[latitude]cos[azimuth], 

Q sin[latitude], 

Q cosjTatitude] sin[azimuth]) 
where Q = 7.29x 1 0‘^ rad/s (rotation of the earth), 

R = 6370320 m (radius of the earth), 

go = 9.80665[l-0.0026373 cos(2xlatitude) +0.0000059 [cos(2x latitude)]^]. 
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The axis system used is as for the point mass model with x being along the line of fire. 

As may be seen from the above equations, the aerodynamic coefficient input 
required is extensive and accuracy of these aerodynamic coefficients is crucial for 
reliable estimates of range. However, for most artillery shells, the reliability of estimated 
values of aerodynamic coefficients is not high enough to inspire confidence m resulting 
range estimates from this model. In spite of this limitation, this model has the capability 
of providing reasonable accoimts of end point data. It can estimate wind corrections and 
can accept non-linear aerodynamic inputs if required. However, because the equation for 
calculating yaw of repose assumes quasi-linear aerodynamics, the use of non-linear 
aerodynamics in equations is questionable. 

3.5 The six-degree-of-freedom model 

As a final step, the most sophisticated trajectory model developed was the six- 
degree-of-freedom model having the 3 spatial degrees of freedom, yaw degree of 
freedom in two planes and the spui. In this model there are no assumptions concerning 
linearised aerodynamics or projectile symmetry. However, the indeterminability of many 
of the initial conditions and aerodynamic coefficients which are required as input 
frequently results in the model not giving significantly better end results than the 
modified point mass model. Thus its usage might not be justified for routine fire control 
work. It is nevertheless a powerful tool for the ammunition designer. 

The six-degree-of-freedom model equations are given as 


du 

7Cpd^ / 

dt 

8m ' 




8m 






(Drag) 

(Lift) 
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where 



^ (Magnus Forces) 



) (Damping Forces) 


go „ r+2(oixu) 
r'^ 

(Gravity and Coriolis) 

(3a) 

• Cl (h-x)vx 

16Ix ’ 

(Spin Damping) 


+ jf. CMav(vxx) 

Oly 

(Overturning Moment) 



(Magnus Moment) 



h X x)] (Damping Moment) 

(3b) 

‘C 

a = cos'~^[(v-x)/v] 




and lyh is the angular momentum vector. 


3.6 Range Tables for B-Shell 

As mentioned earlier, due to non availability of real firing data, the range tables 
for the 155mm Bofors shell HE 77B (here after referred to as B-shell for convenience) 
are used for the present study. These tables and the explanatory notes (given in Appendix 
A) were supplied by ARDE, Pune. The range table lists range obtainable for various 
firing table elevations (hereafter referred to as firiug angle, 0) under standard calm 
atmospheric conditions and for nominal values of weight and muzzle velocity. Also time 
of flight and correction to bearing drift is provided. The corrections to these listed values 
of range for each value of 6 due to following variations at the time of firing are also 
provided: 

1) Variations in ambient atmospheric conditions (temperature, density) 
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2) Head/tailwind 

3) Crosswind 

4) Change in weight and muzzle velocity from the corresponding nominal value. 

The range tables have been prepared using the modified point mass trajectory 

model. As pointed out earlier, this model suffers from a few limitations, mainly because 
of the not-so-reliable accuracy of aerodynamic coefficients that need to be used to solve 
the equations of motion. The complete six-degree-of-freedom model would also suffer 
for the same reasons. Even if accurate aerodynamic coefficients were available, there are 
few other uncertainties that would affect the range estimates in real life and cause 
dispersion of shells. A few such imcertainties are briefly outlined below: 

1. Slight differences between nominally similar projectiles 

Any differences in mass, shape, surfece finish, etc., will cause changes in the 
trajectory of the projectile to some extent. Projectiles might suffer from two types of 
asymmetries while passing through the manufacturing stage: Configurational 
asymmetries and inertial asymmetries. The inertial asymmetries may be of mass 
imbalances wherein center of mass is not on the geometrical axis of symmetry and 
dynamic unbalance where the principal axis of inertia is not coUinear with this axis. 
These asymmetries can lead to large dispersions, and it is by keeping manufiicturing 
tolerances to minimum, and by imparting some amount of spin that one can reduce 
the dispersion due to this effect. However, spin rate is to be carefully chosen to avoid 
the spin-yaw resonance zone. 
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2. Meteorological conditions 

It is noted that the solution of equations of motion uses a single value for head 
or tail wind and crosswind. Furthermore, the corrections to the range (tabulated for 
standard atmospheric conditions) for variation in temperature and density are 
incorporated by using a single value of temperatures and density at the location of 
firing. It is realized that the shell will be passing through different wind and 
atmospheric conditions of temperature and density as its altitude varies during its in- 
flight trajectory. However, for the purpose of applying corrections, a weighted mean 
value of wind velocity, temperature and density is used to account for varying 
conditions prevailing at different altitudes. 

3. Variable launch conditions: Variation in muzzle velocity, gun-jump and throw- 
off, initial yaw, etc. 

At the time of shell leaving the gun barrel, the mitial conditions experienced by 
the shell are not identical for all the shells. In particular, the initial velocity of shell 

(muzzle velocity) depends upon the charge concentration at the time of jSring. 

For the spun projectiles crosswind has an effect in the vertical plane due to the 

time taken to adapt to the air-relative zero yaw position. This occurs whenever the 
projectile emerges into a region of significantly different wind-relative zero yaw attitude, 
but most commonly has its effect on a projectile exiting from the latmcher. Unlike the 
downward carry effect, which is usually roughly quadratic in range and depends on the 
wind at all points down-range, this effect is called aerodynamic jump, which is linear in 
range and dependent on wind at one point. There is a similar term due to the transverse 
component of a head wind. The total effect is called windage jump. 
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The initial yaw and yawing rate of a projectile will obviously determine to a 
large extent what yaw levels are present in the early stage of the trajectory. Furthermore, 
if the damping is poor, they may lead to consistently high yaw and thus cause dispersion 
in range due to yaw drag, especially with indirect fire munitions. 

3.7 The Neural Models 

The limitations of the mathematical models so far used for predicting range of 
artillery shells motivated us to look at an alternative approach to modelling. The feed 
forward neural network provides one such potential way of modelling. The neural model 
needs identification of suitable input-output variables to map the relationship existing 
between them. The functional relation is obtained by training with measured input-output 
variables and then the trained network is used to predict the output for set of inputs not 
seen by network during the training phase. The data firom the range table provided by 
ARDE, Pune for B-shell is used for the present work. 

For the purpose of developing a neural model, the data Jfrom the range table was 
randomly selected to form sets having chosen number of data points, say, 100 or 50 or 25 
or 15. One of these sets is used for training the network. Then two sets of data points, 
each of 30 are taken for validation sets, which are different firom the training set. These 
two sets are used to verify the acceptability of the trained network. As mentioned in 
Appendix B, the acceptability of the network architecture is decided by comparing the 
MSE during training phase and validation phase. The thumb rule applied is that MSE for 
the two validation sets should be not greater than about two times the MSE for the 
training set. Of course, the MSE permitted on the training set is prescribed and kept 
below it, while choosing the architecture of the network, i.e., while choosing the learning 
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rate, the momentum rate, the number of iteration, the number of neuron, etc. Finally, one 
set of data, not used for the above two stages of training or validation is used to predict 
the required output. This output is compared with the known output corresponding to the 
same inputs, and thus shows acceptability of the neural model in predicting the required 
output variables. 

From application point of view, the following three types of modelling 
problems were taken up for study. The neural network modelling for these three cases are 
presented m details in the next chapter along with the results obtained via each of these 
models. However, a brief outline of the three models is given below. 

Model 1 

This model deals with the direct problem of modelling Range, time of flight, and 
drift experienced by the B-Shell. These output variables of the neural model would 
depend on the firing angle 6, air temperature and density, head/tail wind, crosswind, 
muzzle velocity, weight of shell and the initial conditions. 

Model 2 

This model deals with the inverse problem of predicting the firing angle 0 
required to achieve desired range, given the air temperature and density, head/tail wind, 
crosswind, muzzle velocity, weight of the shell. It also predicts time of flight and drift. 
This is what a user (soldier) would require to know in most of the real life applications. 

Model 3 

This model was developed to predict the range that would be obtainable under 
standard conditions, given the measured data of range obtained under varying 
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atmospheric conditions like temperature, density, head/tail wind, crosswind, muzzle 
velocity, and shell weight. Such a prediction capability would be useful to compare 
performance of diflFerent shells under identical (standard) conditions. It is not practical to 
obtain data for different shell or for the same shell on different days under similar 
conditions, and thus compare the relative performance. If data collected under different 
ambient conditions can be used for training the network, and then be able to predict the 
range that would result imder chosen standard ambient conditions, one could then 
compare results for the different shells or the same shell fired at different times or days. 

As a subset of this model, attempt was also made to estimate sensitivity 
coefficients (partial derivatives) that would show how the range is affected by variations 
in each of input variables. These sensitivity coefficients were in turn used to obtain 
standard range from the range data corresponding to ambient conditions of temperature, 
density, wind, etc. 

The results and discussions for all the above three models, along with details of 
modelling are presented in the next chapter. 



CHAPTER 4 

NEURAL MODELS, RESULTS AND DISCUSSION 

In this chapter, we first analyse the data of B-shell supplied by ARDE, Pune. 
Next, details of all the three models mentioned in the previous chapter are presented. 
Finally, the results for all the neural models are discussed. 

4.1 Analysis of B-shell data 

The data sets, explanatory notes and geometric, mass and moment of inertia 
characteristics of B-shell supplied by ARDE, Pune are given in Appendix A. The data 
tables (Table F (i) and Table F (ii)) contain range data for firing table elevation (Firing 
angle) from 0 mil to 1298.4 mils. These firing angles correspond to variation from 0°to 
70° (360° = 6400 mils). Table F (i) has the basic data for range (X, col 1) along with the 
time of flight (t, col 7), correction to bearing for drift due to spin (AcAd, col 8) and 
correction to bearing for 1 knot crosswind (AcAb, col 9) obtainable for varying firing 
angles (Ae, col 2) at sea level under standard atmospheric conditions. Tables F (ii) gives 
corrections to range for the non-standard conditions. Corrections to be applied to 
calculate the range at various firing angles for variation from the nominal values of 
muzzle velocity, projectile weight, air temperature and density, and for mean weighted 
head/tail wind are tabulated. In Table F (ii), column 2 gives correction for muzzle 
velocity being Im/s higher or lower than nominal value, column 3 gives correction for 
Iknot head/tail wind, column 4 and 5 give, respectively, corrections for 1% higher or 
lower value of temperature and density as compared to those for standard atmosphere, 
column 6 gives correction for weight of shell being 1 unit (1 unit = 0.45 kg) higher or 
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lower than its nominal value. The nominal values of the B-shell and its geometry are 
given in fig. A1 of Appendix A. Explanatory notes on the terminology used in Range 
table is given following the range table in Appendix A. 

The range table F (i) and F (ii) were used to prepare input-output data sets for the 
neural models. The range R for non-standard conditions was calculated as follows; 

R = X + Cl (AV) + C 2 (Wx) + C 3 (% BAT) + C 4 (% BAD) + C 5 (APW) (4.1) 
where X = Range at standard condition for a particular firing table elevation ( firing 
angle, 0). 

Cl = Correction to range for muzzle velocity being higher or lower by Im/s than 
the nominal velocity, given in column 2 of the table F(ii) of appendix A., 
sign of Cl is negative if velocity lower and positive if velocity is higher than 
the nominal value. 

AV = Muzzle velocity - nominal value of 818 m/s. 

C 2 = Correction to range for 1 knot head/tail wind, shown m column 3 of table 
F(ii) of Appendix A, sign is negative for head wind and positive for tail 
wind. 

Wx = Weighted mean value of head/tail wind. 

C3 = Correction to range for temperature being 1% higher or lower than value in 
standard atmosphere, given in column 4 of table F(ii) of appendix A, sign is 
negative for temperature being higher and positive for it being lower than 
standard value. 

% BAT = (Actual BAT - BAT in Standard atmosphere) x 100/Actual BAT 

C 4 = Correction to range for density being 1% higher or lower than value in 
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standard atmosphere, given in column 5 of table F(ii), sign is negative for 
density being higher and positive for it being lower than standard value. 

% BAD = (Actual BAD - BAD in Standard atmosphere) x 100/ Actual BAD 

C 5 = Correction to range for 1 unit change in projectile weight, given in the 
column 6 of the table F(ii) of Appendix A, sign is negative for weight 
being higher and positive for being lower than the nominal value of 42.6 kg 

AP W = (Actual weight of Shell - 42.6) 

It is noted that the crosswind gives rise to drift in the lateral direction but does not affect 
the range. 

To illustrate the calculation of range R for non-standard conditions, the following 
example is presented; 

Consider the range for firing angle of 67.2 mils. The Table F(i) gives the 
standard range X = 6600 m, the time of flight =10.2 sec, correction to bearing for drift 
AcAd = 2.8 mils and correction to bearing for 1 knot crosswind AcAb = 0.17 mils. 

Let us assume that the muzzle velocity is 819.9 m/s, tail wind is 8.01 knots, 
crosswind is 8.7 knots, projectile weight is 43.1769 kg, air temperature and density are 
297.3 8 IK and 1.2676 kg/rn^ respectively. For the firing angle of 67.2 mils, from table 
F(i) and F(ii), we have 

Cl = -13.0 m/ m/s 

C 2 = -1.7 m/knot 

C3 = 4.4 m / 1% change in BAT 

C 4 = 17.9 m / 1% change in BAD 

C 5 = 13 m / limit change in projectile weight 
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For the chosen values, the following quantities are easily calculated 
AMuzzle velocity = AV =1.9 m/s, 

AProjectile weight = APW = 1 .282 unit, 

% change in BAT = 3.2, 

% change in BAD = 3.479, 

Substituting above value in Eq. (4.1), we get the non-standard range: 

R=6600 - [(-13.0 X 1.9) + (-1.7 x 8.01) + (4.4 x 3.2) + (17.9 x 3.479) + (13 x 1.282)] 

= 6545.3 m. 

Bearing Correction (T^) = Correction to bearing for drift + Correction to bearing for 

crosswind. 

'P = 2.8 + (8.7 X 0.17) = 4.279 mils 

The implication of the bearing correction is that the projectile would drift to the right 
(looking from behind in the direction of firing) by angle To hit the target, the shell is 
fired by aiming to the left by angle 4^ and this would, in turn, ensure that the target is hit 
after the drift that the shell will ejqperience due to spin and crosswinds. 

Following the above procedure, the data given in the range tables for B-shell was 
converted mto a convenient form for its use in neural models to be developed. 
Specifically, the data yielded sets of data, each set having standard range X, non- 
standard range R, time of flight T, drift angle 'P, ambient temperature and density, 
crosswind, head or tail wind, variation from the nominal value of muzzle velocity and 
projectile weight, or percent variation thereof. 
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4.2 ModeUing 

The data sets generated from the range tables of B-shell were used for training 
the neural models. However, random selection of input-output pairs for any of the 
models (briefly described in chapter 3) showed poor training. A close look at the data 
and considering the physics of the projectile motion gave the clue; as is well known for 
idealized pomt mass projectile motion in vacuum, the range increases as the angle of 
projections increases from 0 to 45 degrees and then decreases as angle increases from 45 
to 90 degrees, the maximum range being at 45 degrees. A similar trend was observable 
for the B-shell and maximum range was just above 45 degrees (802 mils). The 
implication of this feet is that same range is attained for two distinct values of angle of 
firing, one below 45 degree and other above 45 degree. Thus the neural network sees 
inherent contradiction in data and hence the difficulty in mapping the data spanned 
across all the values of firing angles. 

To resolve the above difficulty, the data was partitioned into two bins: one 
having data for &ing angle from 0 to 45 degree and the other from 45° to 70° (firing 
tables gives data up to 70° only). The neural models for these two were separately 
developed. It may be pointed out that, such an approach will not create any difficulty in 
real-life applications because the requirements about the desirable firing angle to be less 
than or more than 45 degree would be known - if the shell is to achieve higher altitude 
during its flight to target, the firing angle of greater than 45 degree would be 
recommended, otherwise less than 45 degree would be the preferred option due to 
shorter time of flight. 
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For the purpose of modelling, a set of randomly selected input-output data was 
selected from each bin. Sets of varying number of input-output samples were tried to 
arrive at the minimum number of samples required for adequate training level. It is 
realized that in real life, the number of samples available will be limited due to the cost 
involved in collecting such range data, and hence the need to search for the min im um 
number of data samples to get an acceptable neural model. The numbers of samples tried 
were 100, 50, 25 and 15. Higher number of samples obviously led to better training, but 
even as few as 15 samples also gave satisfactory models. The suitability of the model is 
tested by a test (validation) set of data, typically consisting of 30 input-output samples, 
selected randomly from the range table; if the MSE was only of the order of two times 
the MSE prescribed for the training phase, the neural model was accepted and its 
architecture fixed for the prediction phase. If not, the architecture was varied till it met 
the above conditions on MSE for training and validation phase. 

Once the neural model is validated, it is used for the prediction purpose. For 
prediction, a set of randomly selected data is taken from the range table, and only the 
input variables of this set are treated as known while the output variable is predicted by 
the neural model. Since the output is also known from the range table, it is compared 
with the neural model output to show how well the predicted values compare with the 
known values. Typically, for prediction, 10 samples were randomly selected from the 
range table. 

In case of neural models like Model 1 and Model 2, prediction of more than one 
output was required, e.g., Model 1 is to predict range, bearing correction and time of 
flight. In view of this, the neural models having more than one output were subjected to 
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the following test: for the same set of inputs, is it better to train the network separately 
for one output at a time, or train it on all the outputs at once. It was observed that even 
though the latter option gave reasonable predictions of all outputs, the former option 
always yielded relatively more accurate predictions. The accuracy being of prime 
concern, it is worthwhile to have separate neural models for each of the output variables 
and, therefore, all studies repeated herein are based on the single output predictions. 

Because the above details of the procedure are common to all the models 
discussed next, we shall refrain from repeating it and only point out wherever the 
procedure differs m some significant way from the above mentioned guidelines. Results 
for the various models discussed next are presented in tabular as well as in graphical 
form. 

4.3 Model 1 

In model 1, we wish to develop neural model for predicting range(R), bearing 
correction ('F) and time of flight (T). These would, therefore, form the output variables 
of the neural model. These output variables depend on the jBring angle (0), ballistic air 
temperature (BAT), ballistic air density (BAD), muzzle velocity (V), projectile weight 
(PW), head/taU wind (Wx) and crosswind (Wz). Thus the neural model has R, or 'R, or T 
as the output variable and 0, BAT, BAD, V, PW, Wx, Wz as the input variables. The 
input variables V, BAT and BAD are given in percent variation from the nominal values 
— nominal muzzle velocity being 818 m/s, and nominal BAT and BAD being sea level 
values in standard atmosphere. The APW value was given as the difference between the 
actual projectile mass and the nominal projectile mass (42.6kg). The schematic of such a 
model is shown in fig. 2. 
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INPUT VARIABLES OUTPUT VARIABLE 


Fig. 2 Schematic Representation of FFNN For model 1 


As mentioned earlier, input-output samples were randomly selected from the 
range table for B-shell. Although training was carried out for data sets having 100, 50, 
25 and 15 samples, the results for the most stringent case of data set having 15 samples 
are presented. For data set having higher number of sanqtles (> 15), the training was 
always superior to that for 15 samples. For illustration, MSB for 15 samples and 100 
samples as a function of iteration is shown in fig.3. It is noted that the minimum MSB is 
achieved fester in case of 100 samples (in less than 200 iterations) as compared to 15 
samples (in about 300 iterations). Furthermore, the value of minimum MSB achieved in 
case of 100 samples is about 8.22xl0‘* where as for the case of 15 samples, it is 
7.64x10"^, which is two orders of magnitude higher than for 100 samples. Fig. 3 also 
shows MSB for the two test data (30 samples) for both cases. It is seen that the MSB for 
both test data is of the order of two times the corresponding MSB values for 15 and 1 00 
samples. 
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Notwithstanding the superior training for the case of 100 samples, for the case of 15 
samples, in absolute terms, the training level achieved was satisfactory, as also was the 
validation test. Hence, all the results presented are only for the case of 1 5 samples. 

For randomly selected 15 samples of 0, BAT, BAD, V, PW, HT and Wz, the 
corresponding R, 'F and T were calculated as explained in 4.1. In addition, two other sets 
of 30 samples each were similarly generated for use as test data. The network tuning 
parameters were varied till acceptable network training was achieved as per criteria set 
out in Appendix B. Now a set of 10 samples, randomly taken from the range table, is 
used to predict R, 'F and T. The results are given m Table 1 and graphically compared 
with the actual values in Fig. 4. As may be seen, the predicted values compare well with 
the actual values. 

4.4 Model 2 

This model is the one, which would be of more use m real-life applications. A 
soldier would want to know the firing angle (0) he should use, the bearing correction ('F) 
he should apply and the time of flight (T) for shell to reach the target. The time of flight 
is required to set the fuse such that the shell explodes after a fixed time - the time being 
chosen to be just prior to shell reaching the target and thereby ensuring that shell 
explodes in air, just before hitting the ground. The information made available prior to 
firing of the shell is the ambient atmospheric conditions: head/tail wind (Wx), crosswind 
(Wz), muzzle velocity (V), projectile weight (PW), atmospheric temperature (BAT) and 
density (BAD). Thus, for the desired range under the prevailmg atmospheric conditions, 
we wish to develop a model that would predict the firing angle, bearing correction and 
time of flight. The schematic of such a model is shown in fig. 5. 
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Table 1. Model 1 : Comparison of Actual and Predicted Range, Bearing 
correction (T^), Time of flight (T) for the given value of 0 


6 in 

Range in m 

^ in Mils 

T in sec 

Mils 

Actual 

Predicted 

Actual 

Predicted 

Actual 

Predicted 

7.7 

1000.8 

1000 

0.097 

0.098 

1.301 

1.3 

11.8 

1498 

1498.5 

0.579 

0.573 

1.898 

1.89 

25.3 

3019.3 

3020 

- 0.071 

- 0.07 

4.119 

4.12 

35.5 

3981.6 

3983 

2.073 

2.075 

5.563 

5.561 

40.9 

4447.8 

4446 

3.557 

3.56 

6.348 

6.34 

52.7 

5424.6 

5423 

3.604 

3.6 

8.049 

8.05 


6406.9 

6406 

5.545 

5.5 

gjjgg 

9.82 


7053.1 

7051 

1.843 

1.845 

11.106 

11.1 

88.4 

8045.9 

8046.5 

■nm 

muH 


13.3 

96.9 

8195.3 

8194.6 

8.991 

9 

13.6 

13.61 
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Fig.5 Schematic Representation of FFNN For model 2 
Before attempting the above model, it was decided to first study simpler 
modelling problem. The number of inputs to the network were restricted to only two: the 
desired range and one from BAT, BAD, AV, APW, Wx, Wz. Fig.6 shows schematic of 
six such submodels studied. It may be noted that, as in Model 1, the BAT and BAD are 
given in percent change from nominal values, and the AV, APW are given in terms of 
difference from the corresponding nominal values. 

The results for Model 2 are shown in Table 2 and Fig.7. The predicted 6, SF, T 
compare well with the actual values, the prediction being marginally better for the lower 
value of 0. One explanation for such observed behaviour is the fact that at lower end of 
0, the range obtainable is more sensitive to variation in 0; typically, increase of 1 mil in 
0 results in range increasing by 134 m for a nominal value of 0 = 0.7 mil, but at nominal 
value of 0 = 802.3 mils, 1 mil increase in 0 would result in range increasing by only by 
3.2m. Notwithstanding this sensitivity of range on nominal value of 0, the predicted 
values by Model 2 are quite acceptable for the whole range of 0 values. 
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Model 2a 



0, or 'i', or T 


Model 2b 


R~ 


% change in BAD 



-► 0,or'F,orT 


Model 2c 



0, or or T 


Model 2d 



>.0, or W, or T 


Model 2e 



>.0, or 'R, or T 


Model 2f 



>0, orT, orT 


Fig.6 Schematic Representation of FFNN for sub models of Model 2 
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Table 2. Model 2: Comparison of Actual and Predicted Firing angle (0), 
Bearing correction for Drift (TO and Time of flight (T) for 
ambient atmospheric conditions, shell weight and muzzle 
velocity considerations. 


Range 

0 in Mils 

Tin Mils 

T in sec 

in m 

Actual 

Predicted 

Actual 

Predicted 

Actual 

Predicted 

20100 

509.10 

508.85 

-0.08 

0.10 

55.50 

57.20 

22800 

738.60 

744.50 

38.60 

37.75 

74.90 

74.80 

23000 

774.60 

773.00 

36.95 

37.10 


76.76 

14000 

233.50 

234.00 

4.10 

4.88 


31.20 

10000 

125.50 

124.40 

10.15 

10.50 

18.00 

17.61 

15100 

272.50 

273.30 

15.04 

15.50 

34.20 

34.29 

12600 

189.70 

191.00 

17.85 

18.80 

25.70 

25.53 

4100 

36.50 

36.40 

0.0 

0.10 

5.80 

5.60 

1100 

8.50 

8.70 

-0.10 

-0.23 

1.40 

1.31 

200 

1.50 

1.62 

0.15 

0.20 

0.20 

- 

0.32 
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Actual beanng corrn 
Pred bearing corrn 


Actual Time 
Pred Time 


Range , mt 


1 1.5 

Range , mt 


1 1-5 

Range , mt 


Fin 7 PomDarison of Actual and Pradictad Thata, Baarino corraction 3nd Tiina of 
Fig. 7 compan \/ dw rat bad. W and W for Model 2. 






The results for the submodels Model 2a to Model 2f are discussed next. It was 
realized that all the submodels could be viewed as a special case of model 2. It would be, 
therefore, of interest to see how the results from the submodels compare with these from 
Model 2 by keeping R and only one of the input variables, corresponding to the 
submodel, nonzero and the rest of the input variables to zero. For example, for submodel 
Model 2a, the Model 2 would have R and % change in BAT as nonzero inputs, and the 
rest of these, i.e., % change in BAD, AV, APW, Wx, and Wz equals to zero during the 
prediction phase. The results so obtainable via submodels and via Model 2 are compared 
in Table 3-8 and Figs. 8-13. The results as a special case of Model 2 corresponding to 
the submodels are given under the heading ’Special - Model 2'. It is noted that, as 
expected, the results via the submodels are uniformly better than those obtained as 
special case of Model 2. 

4.5 Model 3 

It is of interest to develop a model that can yield standard range (X) obtainable 
under standard conditions - temperatures and density corresponding to sea level values 
in standard atmosphere, no head/tail wind, shell having standard muzzle velocity and 
weight. The data assumed to be available for such modelling is the measured range (R) 
for existing ambient temperature, density, head/tail wind, muzzle velocity and shell 
weight. Two distinct approaches are used and the models, named model 3 A and model 
3B, have been developed for the purpose of predicting standard range. 

Model 3A 

The standard range differs from the non-standard range due to the effect of i) 
variations in atmospheric temperature and density, ii) amount of head/tail wind, iii) 
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Table 3. Model 2a: Comparison of Actual And Predicted Angle of Firing (0), Bearing Correction for 
Drift (¥) and Time of Flight (T) for Ballistic air temperature (BAT) Consideration. 
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Time of flight , sec Bearing correction , mils 





Fig. 8 Comparison of Actual and Predicted Theta, Bearing correction, and time of 

flight for varying Range and BAT. 
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Table 4. Model 2b: Comparison of Actual And Predicted Angle of Firing (0), Bearing Correction for 
Drift ('T) and Time of Flight (T) for Ballistic air density Consideration. 
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Table 5. Model 2c: Comparison of Actual And Predicted Angie of Firing (0), Bearing Correction for 
Drift ('P) and Time of Flight (T) for Velocity Consideration. 
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Bearing correction , mils 





Fia. 10 Comparision of Actual and Predicted Theta, Bearing correction, 
and Time of flight for varying Range and V. 
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Table 6. Model 2d: Comparison of Actual And Predicted Angle of Firing (0), Bearing Correction for 
Drift Q¥) and Time of Flight (T) for Projectile weight Consideration. 
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Time of flight , sec Bearing correction , mils 
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X 10'' 



Fig. 11 Comparison of Actual and Predicted Theta, Bearing correction, and time of 

flight for varying Range and PW. 
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Table 7. Model 2e: Comparison of Actual and Predicted Angle of firing (0), Bearing Correction for 
drift; (T^) and time of flight (T) for head/tail wind consideration. 
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Table 8. Model 2f: Comparison of Actual and Predicted Angle of Firing (0), Bearing Correction for 
Drift ('P) and Time of Flight (T) for crosswind consideration. 
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Bearing correction , mils 









variations in muzzle velocity and shell weight from the nominal values. If one could 
evaluate sensitivity coefficients that estimate how much the range varies per unit change 
in each of these, then one can calculate the standard from the non-standard range as 
follows: 


X = R-f®VT-(®lAp, 

I ST j I ap J 


^Iav- 

svj Ispwj 


where X is the standard range, R the non-standard range. 






are, respectively, the sensitivity coefficients with respect to 


temperature, density, head/tail wind, muzzle velocity and shell weight. AT = T - To, Ap = 
p - Po, AV = V - Vo, APW = PW - PWo, where T, p, V, PW are the values of temperature, 
density, muzzle velocity and shell weight at time of firing, and To, Po are the sea level 
standard values, Vo and PWo are the nominal value of muzzle velocity and shell weight. 
The Wx is the value of head or tail wind (Wx is positive for tail wind and negative for 
head wind). To estimate sensitivity coefficients, five separate neural models (shown in 
fig. 14) were developed. The methodology used is based on the Delta method proposed in 
ref. [8,13] for estimating aircraft stability and control derivatives from the flight data. The 


derivatives (sensitivity coefficients 




etc.) can also be estimated by using 


the Delta method. Here the FFNN is trained to map one of the input variables from T, p, 
Wx, V, PW to the output variable range(R). The network input is now given a small 
perturbation in both increasing and decreasing directions. Such a modified input file is 
now presented to the trained network to predict the Range, R at its output node. The 


50 



difference in the predicted value of the range from the value predicted for the original 
value of the input is attributed to the perturbation of the chosen input variable. The 
difference so calculated in the value of range divided by the perturbation value 
introduced in the input variable yields the corresponding sensitivity coefficients. For 


example, to estimate 




the temperature T is varied to T + AT and T - AT and the 


corresponding predicted range R'*' and R' are recorded. Then 
f np \ R"^ — R~ 

— = , Por the purpose of illustration, non-standard range corresponding to 

ar j 2AT 

standard range of 20000 m was calculated for various combinations of T, p, Wx, V and 
PW. A set of 15 such data was used to train the five neural networks having non-standard 
range as the output and one out of T, p, Wx, V, PW as the input. The trained network was 
used to estimate the five required sensitivity coefficients; one from each of the networks 
shown in fig. 14. Finally, the standard range is calculated by using Eq. 4.2. A conqiarison 
of predicted and actual standard range is given in Table 9. Table 9 lists the non-standard 
range along with the sensitivity coefficients. Table 9 shows that the standard range 
estimates via sensitivity coefficients compare well with the actual standard range. 
However, this approach is a round about way of predicting standard range, and, therefore, 
a more direct method would be highly desirable. This search for a direct method led to 

the following model. 

Model 3B 

Model 3B has two versions, named Model 3B-1 and Model 3B-2 that were 
studied. Model 3B-1 is identical to Model 1, but it is used in a special way for estimating 
standard range. Like Model 1, output of the network is range R, which is mapped to 9, 
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Pig.l4 Schematic of Neural Models for Estimating Sensitivity Coefficients 
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Table 9. Model 3 A: Comparison between actual and calculated range at standard atmospheric 

conditions. 
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AV, %BAT, %BAD, APW and Head/tail wind (Wx). However, mapping of the inputs to 
bearing correction ^ and time of flight T is not required as done for Model 1. The 
network is trained on randomly selected 40 data taken from the range table. Now, the 
trained network is used to predict standard range for each of the firing angle 0 by setting 
AV = APW = %BAT = %BAD = Wx = 0. The standard range for various value of 0 so 
obtained is given in table 10. 

Model 3B-2 is similar to the Model 3B-1 except that instead of AV and APW, 
percent change from the nominal value of the muzzle velocity and the shell weight were 
used as the input variables. The trained network was again used to predict standard range 
for various firing angles by putting all variables except 0 to zero. The results are shown in 
Table 10. 

In Table 10 and Fig.15, the predicted standard R via Model 3B-1 and Model 3B-2 
are compared with that from the range table. Though both the models yield standard R 
which compares well with that given by the range table, it is observed that the Model 3B- 
2 gives relatively better results compared to Model 3B-1. Also, it may be noted that 
predictions for standard range are relatively better for lower values of 0, i.e., there is 

slight decrease in prediction accuracy of standard R as 0 increases. 

It is emphasized that both options of Model 3B are more suitable and 
straightforward to use for real life applications. In contrast, model 3A requires multiple 
networks to be used and calculation of the standard range using results from these 

networks. 
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Table 10. Comparison of Actual Range at standard condition for 
Different types of model. 


0 in 
degree 

R-non 
standard 
in m 

X-s-bandard 
in m 

X- Model 
3B-1 in m 

X- Model 
3B-2 in m 

25.41 

18922.00 

19100,00 

19107.00 

19101.00 

25.72 

19315.00 

19200.00 

19205.00 

19199.00 

26.03 

18921 . 00 

19300.00 

19304.00 

19298-00 

26.34 

19439.00 

19400.00 

19404.00 

19399.00 

26. 66 

19688.00 

19500.00 

19503.00 

19498.00 

26. 98 

19097.00 

19600.00 

19602.00 

19598.00 

27.30 

20226.00 

19700.00 

19701.00 

19696.00 

27.62 

19900.00 

19800.00 

19800.00 

19796.00 

27.96 

19442.00 

19900.00 

19900.00 

19896.00 

28.29 

19996.00 

20000.00 

20001.00 

19997.00 

28.64 

19861.00 

20100.00 

20101.00 

20097.00 

28.98 

19991.00 

20200.00 

20201.00 

20197.00 

29.33 

20545.00 

20300.00 

20302.00 

20299.00 

29.69 

20667.00 

20400.00 

20402.00 

20399.00 

30.05 

19922.00 

20500.00 

20502.00 

20499.00 

30.42 

20355.00 

20600.00 

20603.00 

20600.00 

30.79 

20748.00 

20700.00 

20703.00 

20700.00 

31.17 

20493.00 

20800.00 

20804.00 

20801.00 

31.56 

21139.00 

20900.00 

20903.00 

20900.00 

31.95 

20933.00 

21000.00 

21003.00 

21001.00 

' 32.36 

21505.00 

21100.00 

21104.00 

21102.00 

32.77 

20819.00 

21200.00 

21204.00 

21202.00 

33.18 

21467.00 

21300.00 

21303.00 

21301.00 

33.62 

21531.00 

21400.00 

21403.00 

21402.00 

34.05 

21678.00 

21500.00 

21503.00 

21502.00 

34.50 

21394.00 

21600.00 

21602,00 

21602.00 

34.97 

21627.00 

21700.00 

21701.00 

21701.00 

35.44 

21738.00 

21800.00 

21800.00 

21800.00 

35.93 

21625.00 

21900.00 

21898.00 

21898.00 

36.44 

22254.00 

22000.00 

21997.00 

21998.00 

36.96 

22141.00 

22100.00 

22095.00 

22096.00 

37 . 50 

22003.00 

22200.00 

22192.00 

22193.00 

38.07 

22370.00 

22300.00 

22291.00 

22292.00 

38.67 

23179.00 

22400.00 

22390.00 

22391.00 

39.30 

22571.00 

22500.00 

22488.00 

22490.00 

39.99 

22071.00 

22600.00 

22590.00 

22591.00 

40.73 

22452.00 

22700.00 

22692.00 

22693.00 

41.55 

23353.00 

22800.00 

22796.00 

22797.00 

42.46 

23171.00 

22900.00 

22903.00 

22902.00 

43.59 

22827.00 

23000.00 

23021.00 

23017.00 
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Fig. 15 Comparison of Actual and Predicted Standard Range for different 

types of model 
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CHAPTER 5 


CONCLUSION 

5.1 Conclusions 

The validity of the neural modelling is demonstrated for three different 
applications relevant for artillery shells. A set of available firing data in terms of range, 
firing angle, shell weight, muzzle velocity, existing atmospheric conditions like 
temperature, density, head/tail wind, cross wind, time of flight, bearing correction are 
required to train the network with suitable input-output samples: which of these measured 
variables would form the inputs and the outputs is decided by the purpose of the neural 
model. In particular, for three distinct applications, a neural model has been identified 
and validated for B-Shell data. The results for all the models compare well with the 
known results. The strength of the neural models lies in the fact that once the network is 
trained, it can be used for on-line applications on the field of action - to predict the range 
for chosen angle of firing or to predict the firing angle for desired range, under the 
prevailing atmospheric conditions and for known shell weight and muzzle velocity. 

The conventional approach, i.e., the mathematical models such as in-vacuo 
trajectory model, modified point mass model, sk-degree-of-fi:eedom model require 
knowledge of all the forces and moments acting on the shell. Evaluation of forces and 
moments, in turn, require aerodynamic coefficients as inputs and this fact limits the 
accuracy of predictions because the reliability of available estimates of these coefficients 
is not always high. In contrast, the proposed neural models do not require any 
mathematical model or its solution. This implies that the neural models do not require 
estimates of aerodynamic coefficients. Furthermore, if the neural network is trained on 
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the real data, it will automatically account for the initial conditions in an implicit way. 
But the mathematical models have no provisions to account for initial conditions such as 
gun jump and throw-off. 

5.2 Suggestions for Future work 

1) The training of the neural network depends on the architecture of the FFNNs, the 
activation function used by the neurons and the values selected for the tuning 
parameters. From the neural network point of view, lot of scope exists to explore and 
experiment with new evolving schemes and methodologies for improving the input- 
output mapping. 

2) The present work has used the firing tables provided by ARDE, Pune. An atten^t 
should be made to seek real firing data and validate the present models for such data. 
Using real data, it will be of interest to predict range for known firing angle, or to 
predict firing angle required for desired range using the mathematical models like the 
modified point mass model or six-degree-of-fireedom model and also from the 
proposed neural Model 1 or Model 2. A comparison of results from mathematical 
models and neural models would show relative reliability of the predicted value for 
real life apphcations. 

3) The data set supplied by ARDE, Pune is valid for weighted value of the head/tail 
wind and cross wind. This is an approximation to varying wind conditions at different 
altitudes through which the shell passes. It may be worthwhile to search for a way to 
account for varying wind conditions. The same comment also applies for finding a 
way to account for varying temperature and density conditions. 
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Table F(i) Table F(ii) 

CHARGE 9 -SHELL 77B CHARGE 9 - SHELL 77B 

Basic Data Corrections to race: for non-standard conditions 
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Explanatory Notes on 155mm BOFORS Range Tables fRT^ 
Basie Data (Table F (i) ) 

y-V t Okl_ T ' ^ _X-£._i 


Column No, 

— — i ^ lU-l 

Description 

1 

Range (X) ; Range in m, being aimed 

2 

Firing Table Elevation (Ag)' Elevation in mils (360 deg = 6400 mils) to be 
set on the gun to achieve range at col. 1 . 

3 

Fuze setting (FS) setting to be done on the fuze so that it 
functions(burst 3 during flight) at the time mentioned at col. 4 ( in case 
time fuze being used) 

4 

Corrn to fuze setting(FS) to change height of burst by dovvn 10 m. 

5 

Effect OQ range, in m, for increase of 1 mil in firing table elevation 

6 

Fork (also known as short bracket) : when a large no of shells are fired 
at the MV and elevation in a series (i.e. in quick succession) they will 
exhibit some dispersion pattern, which (supposedly) follows normal 
statistical distribution. The x(range) and y (cross range) co-ordinates of 
fall of shot are treated separately as independent variables and the 
standard deviation Ox & oy in x & y resp. are obtained 50 % zone in 
range and cross range is defined as the length zone spread equally 
around the mean point of impact (more commonly known as MPl) of a 
series of rounds fired which may contain 50% of the rounds if a large no. 
of rounds are fired. 50% Zone = 1.349 ox (or oy as applicable) Half the 
50% zone is known as probable error, i.e 50% zone consists of one 
probable error on either side of the MPl. Fork : the change in elevation 
required to produce a range change equivalent to 4 x range probable, 
errors. This parameter is mainly used for range safety purposes. 

7 

Time of flight : the time taken by the projectile to travel from the 
moments it leaves the gun barrel to the point of impact. 

8 

Correction to Bearing for drift' The angular correction to be applied to 
the line of fire called bearing (bearing is the angle measured clockwise 
from local North) to compensate for the deviation of the shell from the 
line of fire towards right due to its right hand spin. Most of the artillery 
shell are gyro, stabilized hence have an eq'"^ total angle of attack (called 
yaw of repose in ballistics) towards right and upwards wrt the velocity 
vector resulting in the deviation of projectile to right of line of fire. 

9 

Correction to bearing for 1 knot cross wind ; corrn to be provided to the 
bearing of line of fire to nullify the deviation of the shell due to constant 
cross wind of 1 Knot magnitude through out the trajectory at all heights. 
The correction is - ve if the crosswind is blowing from Lt to Rt i.e. + ve 
cross wind In practice the wind will never be constant, at all heights. A 
weighted mean value (wrt heights) is obtained called as equivalent 
constant wind (ECW) A wind of quantity ECW Vvill produce the same 
effect as that of the actual-prevailing wind. Hence tor the same prevailing 
wind conditions ECW values will be different for different ranges due to 
the difference in the vertex heights. Same thing is applicable in case of 
range i e. head or tail wind The measured wind data at different is first 
resolved in to range and cross wind components at all heights. 
Equivalent constant winds for range and cross components are 
calculated and used separately. 

10 

Corrn to "corrector" to change height of burst by down 10 m. 

t 


j 
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Correction to range for non-standard conditions (Table F (ii) ) ; 


Column No. 

Description 


I 

2 & 3 

I 

t 

f^nngo (X) : Rat\go in m, being aimed 

Correction to range in m for a change in muzzie velocity (MV) by 1 m/s 
f the MV increases then the range also increases, hence the correction 
.0 range for increase in MV is -ve 

4&5 

Correction to range in m for a range wind of magnitude 1 knot 

A tail wind is taken as +ve & head wind -ve 

f it is tail wind then the range increases, hence the correction to range ^is 
-ve 

6 & 7 

Correction to range in m for a change in ballistic air temp by 1% of the 
standard . 

The RT is compiled for some standard meteorological conditions, which 
defines the air temp , pressure and density at various altitudes. In case of 
155mm RT ICAO standard atmosphere is assumed 

Air is assumed to be dry in ICAO but in practice it may have some 
moisture. Also the temp and press may not be the same as the 
standard 

Ballistic air temp. (BAT) ; BAT is the temp, of the dry air having,an equal 
density as that of moist air. The correction tor BAT only takes care of the 
Changes in range on account of the changes in Mach No. due to variation 
in temp, (recall that acuastic velocity = V(T RT) ) and not the variation in 
density due to changes in temp. This is also termed as elasticity effects. 

8 & 9 

Correction to range in m for a change in ballistic density (p) or 
change in drag coefficient (Colby 1 %. If p or Co increases then the 
range decreases, hence the correction to range is +ve. 

■ Density at different altitudes will be different, however a single weighted 
mean value is obtained ( in the same manner in which wind is treated). 

10 & 11 

Projectile weight : Correction to range for a projectile whose weight 
differs from the standard weight by 1 unit (□ mpj : D marking is 
done on the shells after manufacture for there identification in 
terms of weight class. One such square means that the shell weight 
differs from the std. by 1 unit ). The method of marking is different 
for different equipments,. 

The effect of change in shell weight is two fold: first the MV reduces 
resulting in the reduction in range and second the carrying capacity 
increases (Drag / Mass decreases i.e. less deceleration than the normal) 
resulting in increase in range. The net effect could be either way. 

For shell HE 77B 1 (□ mp!)= 0.450 kg . 1 


Explanation about “Ammunition Data(cont.)” on page (x) 

Shell weight influence:- 

The relationship given there viz. A Mv = A2 x A 1 00 g shell weight, 

Gives the change(decrease) in muzzle velocity (Mv) due to 1-00 gm change(increase) 
in the shell weight for various charges (M3A1 3G etc are designated propellant 
charaes with fixed qty of propellant. Each charge corresponds to a certa n Mv. It can 
be seen that the shell influence is increasing as the Mv increases (from top to 


bottom) 
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Nose 

Total Length = 825 mm 

CG Location = 533 mm from Nose-tip 

Moment of Inertia- 

Axial, Ixx = 0.146 kg-m^ 

Transverse, lyy “ 1.709 kg-m^ (through an axis passing through CG) 

Twist of Rifling(n) = 1 in 20 calibres. the shell will complete 

one full rotation (360°), in side the gun barrel by 
the time it does a translation or linear travel of 20 
calibres i.e. 20 * 0.155 = 3.1 m. 

Thus launch spin = Mv / (n.d) = 818 / (20 * 0.155) = 263.87 rev/s = 15832 rpm 

Fig A1 . Schematic of 155 mm Bofors Shell HE 77B 
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Appendix B 

FEED FORWARD NEURAL NETWORK 

The back propagation network consists of one input layer, one output layer and 

one or more hidden layers. There is no theoretical limit on the number of hidden layers 
but typically there is just one or two. Some work has been done [14], which indicates that 
minimums of four layers (three hidden layers and one output layer) are required to solve 
problems of any complexities. Each layer is fully connected to the succeedmg layer 
(standard connection). 

There are as many neurons in the input layer as there are inputs, and likewise with 
the output layer. The number of layers and the number of neurons in the hidden layer(s) 
must be determined by trial and error. There is no quantifiable best answer to the layout 
of the network for any particular application. There are only general rales picked up over 
time and followed by most researchers and engineers applying various architectures to 
their problems. 

Rule 1. As the complexities in the relationship between the input data and the desired 
output increases, the number of processing elements in the hidden layer should increase. 

Rule 2. If the process being modeled is separable into multiple stages, then additional 
hidden layers may be required. If the process is not separable into stages, then additional 
layers may simply enable memori 2 ation and not a true general solution. 

Rule 3. The amount of available training data sets an upper bound for the number of 
processing elements in the hidden layer. To calculate this upper bound, use the number of 
input-output pair examples in the training set and divide that number by the total number 
of input and output processing elements in the network. Then divide the result again by a 
scaling fector between five and ten. Larger scaling factors are used for noisy data. 
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Extremely noisy data may require a factor of twenty or even fifty. Very clean input data 
with an exact relationship to the output might allow the factor to be dropped to around 
two. 

Cross-Validation and Overtraining 

One approach to avoid over-training of the network is to estimate the 
generalization ability during training and stop when it begins to decrease. The essence of 
back-propagation learning is to encode an input-output relation, presented by a set of 
data, with a multiplayer perceptron well trained in the sense that it learns enough about 
the past to generalize to the future. The simplest method is to randomly partition the data 
set into a training set and a test (validation) set. From the tra ining set, a validation subset, 
which are typically 10 to 20 percent of the training set is set aside. The motivation here to 
validate the model on a data set different from the training set that is used for selecting 
the architecture of the network. The training set is used to modify the weights, the 
validation set is used to estimate the generalization ability. The architecture of the 
network is varied till the training set results in MSE less than the prescribed value e. This 
architecture is now tested on the test data (which can be one or more) and if the MSE is 
of the order of 2e, the architecture is accepted to yield the desired neural model and 
assumed to be capable of predicting required output for iiq)uts not seen earlier by the 
network 

Another way of avoiding over-training is to limit the ability of the network to take 
advantage of spurious correlation in the data. Over fitting is thought to happen when the 
network has more degrees-of-freedom (the number of weights, roughly) than the number 
of the tr aining samp les when there are not enough examples to constrain the network. 


75 



I* n I l!’ » 


Even though it may give exactly right output at the training points, it may be very 
inaccurate at other points. An example is a higher order polynomial fitted through a small 
number of points. 

Sufficient Training Set Size For a VaKd Generalization 

Generalization is influenced by three factors: i) the size and the efficiency of the 
training set, ii) the architecture of the network, and iii) the physical complexities of the 
problem at hand. Clearly, we have no control over the last factor, i.e., the physical 
complexity. We have already discussed the choice of architecture based on training and 
test data. Once the architecture of the network is fixed, then the size of training set can be 
derived as follows. 

Let M denote the total nmnber of hidden layer computation nodes. Let W and N 
be the total number of synaptic weights and the number of random examples used to train 
the network respectively. Let 8 denote the fraction of error permitted on test. Then, 
according to Baum and Haussler, [2] the network will almost certainly provide 
generalization provided the following two conditions are met. 

(a) The fraction of error made on the trainmg set is less than s/2. 

(b) The number of examples(N) used in the training is 


W W 
N>32 — In — 
s s 


where W is the total number of synaptic weights. 

Ignoring the logarithmic factor, taking first order approximation, the number of training 
examples is directly proportional to the number of weights in the network and inversely 

W 

proportional to the accuracy parameter s. Then, N > — . 


76 



